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ABSTRACT 

The blade root response due to ice impact on an engine blade is simulated using the NASA in- 
house code BLASIM. The ice piece is modeled as an equivalent spherical object impacting on the 
leading edge of the blade and has the velocity opposite to that of the aircraft with direction parallel 
to the engine axis. The effect of ice impact is considered to be an impulse load on the blade with its 
amplitude computed based on the momentum transfer principle. The blade response due to the 
impact is carried out by modal superposition using the first three modes. The maximum dynamic 
stresses at the blade root are computed at the quarter cycle of the first natural frequency. A 
combined stress failure function based on modified distortion energy is used to study the spanwise 
bending damage response at the blade root. The damage function reaches maximum value for very 
low ice speeds and increases steeply with increases in engine speed. 


INTRODUCTION 

At high altitudes, when aircraft flies through clouds of super-cooled water droplets, ice formation 
occurs on forward facing structural components. One such component is the engine inlet. With 
time, the ice accretes on the inlet and eventually sheds due to structural vibrations. As a result, 
blocks of ice travelling at the speed of aircraft impacts the engine blades rotating at high RPM. A 
schematic of this phenomenon is shown in Figure 1. This process may cause severe damage at the 
blade root and subsequently may result in catastrophic failure of the engine. 

In this paper, the blade root response due to ice impact on an engine blade is simulated using the 
NASA in-house code BLASIM [1]. The ice piece is modeled as an equivalent spherical object 
impacting on the leading edge of the blade. The ice piece velocity is opposite to that of the aircraft 
with direction parallel to the engine axis. The effect of ice impact is considered to be same as that 
caused by an impulse load on the leading edge of the blade with an equivalent amplitude computed 
based on the momentum transfer principle. The dynamic analysis is carried out by finite element 
method andthe blade response due to this impact is carried out by modal superposition. However, 
experience has shown that the highest root stresses occur at the quarter cycle of the first bending 
frequency [2]. A combined stress failure function based on modified distortion energy [3] is then 
used as a measure of the root damage response due to impact. The details of the linear response 
analysis as implemented in the code are presented. The root damage response capability of the 
code is demonstrated by considering a modified SR-2 unswept propfan blade made. Parametric 
studies are performed and the effects of ice speed, ice size, engine speed and impact location on the 
root damage response are discussed. 
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ANALYSIS 


Ice Impact Model 

The geometry of ice impact on the leading edge of the blade is shown in Figures 2(a,b). The ice 
piece is modeled as an equivalent spherical object impacting on the leading edge of the blade. The 
impact on the blade arises from the momentum of the ice piece. The duration of impact is 
considerably short and for the purpose of obtaining the response at the root, the impact is 
considered as an impulse on a single node along the leading edge of the blade. Also, the impulse 
component normal to the chord, being the primary factor in causing the spanwise bending damage, 
is included in the response computation. The impact velocity direction relative to the blade is a 
function of the aircraft speed and the rotational speed of the blade (Figure 2a). The resulting impact 
depends on the magnitude and direction of the relative velocity and the mass of ice piece. 
Depending on the spacing (i.e., number) of blades and relative velocity, the ice piece that is 
approaching the blade under consideration is sheared off by the adjacent blade as shown in Figure 
2b. As a result, only a portion of ice piece impacts the leading edge of the blade. Effectively, the 
size of the ice piece that finally impacts the blade depends on the blade spacing, blade speed and 
aircraft speed. The shape of the ice is considered to be a sphere only for simplicity and affects the 
determination of sheared part of the ice. Once the initial velocity equivalent to the impulse is 
computed, the remaining analysis does not depend on the ice shape. 

Blade Finite Element Model 

In the code, the blade geometry is input in the form of finite element grid and nodal thicknesses. 
The blade is modeled using three node triangular plate finite elements with six degrees of freedom 
per node. This finite element is similar to NASTRAN (TRIA3) [4,5]. The blade finite element 
mesh consists of 55 nodes and 80 elements. The planform of the blade considered for analysis and 
the corresponding finite element blade mode are shown in Figures 3a and 3b. 

Response Analysis 

The blade root response analysis, as implemented in the code is presented here. The governing 
equations of motion for an undamped discrete system of order n is given by 

[m]{q\ + [k][q] = {f} ( 1 ) 

where [m] is the mass matrix, [/t] is the stiffness matrix, {< 7 } is the vector of n degrees of freedom 
and {/} is the forcing function. 

Consider the first N modes (N<n) and let [u] be the modal matrix of the above undamped system 
with N eigenvectors as its columns, i.e., 


[u] = [{«}U«}Mk}*] (2) 

and {<7) = [M]{T7«} (3) 

then Equation (1) transforms to 

[M]{v} + mv) = M T tf) (4) 

where [ M] and [K] are the generalized mass and stiffness matrices, tj and ij are the generalized 
coorduinate and the associated acceleration vectors. 

The matrices [M] and [K] are diagonal due to orthogonal property of the modal matrix [«]. 
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Premultipling the Equation (4) by [M]' 1 , 

{ij} + [tt> 2 ]{>7) = IW] 1 [«] T 0 r ) (5) 

where <o's are the natural frequencies of the system and [g£].= \M\‘^\K\. 

Since the matrix Equation (5) is uncoupled, it can be written as 

n 

7]r(0+ (#T]r(i) = tJt X U r Ji(t) , T = 1,.JV (6) 

M r i=l 

where u\ is the 1 th component of the r* eigenvector and M r is the generalized mass for the H* 1 
mode. If the forcing function is acting only on one of the degrees of freedom, say i-j, then the 
Equation (6) simplifies to 

ij r (t) + 6Js?7} r (r) = F r (t) , r=\,..N ( 7 ) 


where F r (t) = jjrfj(t) 

Instead of the forcing function, F r (r), if an impulse acts on the r* generalized degree of freedom, 
the set of equations can be solved by applying an initial velocity condition, i.e., 

jjr(t) + CDr 2 T]r(t) = 0, (r = 1 ,.M ) with T] r (0) = F r (8) 


^ l^r ^ 

where F r = and f, is the magnitude of impulse on y* degree of freedom. 
M r J 

The solution of equation (8) can be written as 


U; * 

rir(t) = rf-fj sinco r t, (r = 1 ,..N ) 

M r COr J 


(9) 


Now, using the relation {q} = [u][r}(t)) y the nodal displacement response can be obtained as 

N N 


U; 


qi(t) = X u\T]r{i) = X uljp-f j sinCDrt , (i=l,..n) 

r=l r*= 1 rrlfOjfr 


( 10 ) 


The q's can also be expressed in the form 

{ q } = Tji(t){u} 1 + tl2(0{«} 2 + m(0{«) 3 + — + T] N (t)[u} N (11) 

if{S!' is the element stress vector that correspond to the mode shape {m}\ then total stress response, 
{<7} can be written as 

(o(0) = + *)2(0{<7) 2+ Tf 3 (t)[c} 3 + .... + T} N (t)[a) N (12) 
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The response of the blade due to impact (impulse loading) is expected to be predominantly first 
mode and hence the maximum stress occurs at or near the quarter cycle of the first mode frequency 

[2]. The stress vector at time, t = nl2(Oi is 

N 

= X (13) 

2 ( 0 ] »=i 2 ( 0 \ 


Computation of the Root Damage 

The BLASIM code is capable of analyzing metallic as well as composite blades. In order to account 
for the complex loads and microscopic failure mechanisms, a combined stress function based on 
the modified distortion energy [3] is used as a measure of the blade failure. The general expression 
for a composite ply, g, is given by 

-feiMg)’- 


0110^220 " 12 * 


(14) 


where a - 


{ 


T (for tension) if <7n £ 0 
C (for compression) if C\ i 


< 0 


P = 


T (for tension) if 022 ^ 0 
C (for compression) if <722 < 0 


. (1 + 4v )2 - Vl3)E22 + (1 - V23)Ell 

12/J nP [El 1^22(2 + V 12 + v 13 )(2 + V 21 + V23)] 1/2 

S is the strength, E is modulus of elasticity, and V is Poisson's ratio with subscripts 1 and 2 
referring to the longitudinal and transverse directions of the fiber in a ply. 

It is assumed that y ' 12 ^ = 1.0 and v 23 = i.25v 12 in implementing the above equation in the code as 
these are typical values for a normal fiber composite system. The values of g are computed for all 
the elements (second row) at the blade root and the maximum is presented as a measure of blade 
root damage in this paper. Material is considered failed when the function value is equal to or 
greater than unity. 


RESULTS AND DISCUSSION 

The SR2 model unswept titanium blade [7] is considered to study the effect of ice impact on the 
root damage as shown in Figure 3a. The corresponding discretized finite element model is shown 
in Figure 3b. The setting angle (orientation of the blade chord with respect to the plane of rotation 
at 75% span) is 57° and the number of blades are eight. The titanium properties used in the 

numerical computation are: E = 16.5xl0 6 pri, G - 6.4x1 0 6 pri, V = 0.3, p = 0.00044 Ibfsec^lin 2 , 
S 11 = 7.4X10 4 psi, 522= 7.4xl0 4 pri,>5i2= 4.4x10* psi. The range of engine speed is 500 - 
10000 RPM and the maximum ice radius considered is 0.8". 

The dynamic deformation shape of the blade is expected to be a simple one because of an impulse 
load at a point. Hence, the response computation through modal superposition technique is carried 
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out using first three modes. The midplane of the 3-D blade and the first three mode shapes (first 
bending, second bending and first torsion) are shown in Figures 4(a-d). In order to have a better 
3-D visual effect, nodal data is interpolated and a dense mesh is shown in these Figures. 

Effect of Ice Size and Speed. 

The effects of ice size and ice speed at 3000, 5000 and 8000 RPM on the failure function are 
shown in Figures 5(a-c). The root damage is maximum when the ice size is maximum at all the 
engine speeds as expected. However, for a specified ice size, the damage reaches a maximum at 
relatively low ice speeds and thereafter decreases continuously. The optimum ice speeds that 
results in maximum root damage for the three engine speeds are 20, 30 and 40 knots respectively. 
The corresponding failure function values are 0.026, 0.076 and 1 .22. This behavior is due to the 
fact that only the impulse component normal to the blade chord is considered to cause the 
spanwise bending damage. At a given engine speed, the impact relative velocity on the blade 
increases with increase in ice speed. As a result, the impulse on the blade increases but the impulse 
component normal to the blade chord decreases due to decrease in the impact angle relative to the 

chord, 6, as shown in Figure 6. For larger ice speeds, the value of 6 gradually reduces to zero (see 
also Figure 2a). 

Effect of Engine Speed 

Figures 7 and 8 show the effect of engine speed on the root damage. The graphs shown are for an 
ice size of 0.8" radius and the impact location at 73.5% span along the leading edge. It can be seen 
from Figure 7 that the failure function increases steeply for larger values of engine speed ( 5000 - 
8500 RPM) In the speed range approximately 8500 - 9000 RPM, it decreases steeply. The damage 
function, as expressed by the Equation (15) is a function of combined stresses due to centrifugal 

loading and impulse loading. Furthermore, centrifugal loading is ~ I2 2 and hence contributes to the 
damage more at higher values than that due to impulse. 

The steep decrease can be attributed to the shearing action of the neighboring blade due to its 
proximity. In the absence of shearing action, the ice mass remains constant and impulse increases 
linearly with engine speed as shown by dashed line (Figure 8). However, due to shearing, the 
impacting mass of the ice and consequently the impulse reduces as shown by solid line. The ice, 
as the engine speed increases, start getting sheared off at approximately 5800 RPM by the adjacent 
blade (see Figure 2b). 

Effect of Ice Impact Location 

The impacting location of the ice is varied from approximately 30% to 75% of the blade span to 
study its effect on the root damage. Figure 9 shows that the relationship between the magnitude of 
root damage and the impact location is almost linear and the slope or the line depends upon the 
engine speed. As discussed in the previous section, the engine speed has a strong effect on the root 
response and hence the slope of the line for higher engine speeds is higher. For example, at 8000 
RPM, the failure function has a value of 0.33 at 30% span and increases to 1.00 at 75% span. 


CONCLUSIONS 

The response at the root of an SR2 titanium blade due to ice impact is simulated via modal 
superposition technique. The failure function based on Modified Distortion Energy criteria is used 
as a measure of root damage. The following conclusions can be drawn from the results presented 
here: 
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1 . The blade root damage due to spanwise bending increases with the increase in ice size and is 
maximum for low values of ice speeds. 

2. At engine speed larger compared to ice speed (aircraft speed), ice piece gets sheared off 
resulting in lower impact mass and consequently low impulse on the blade. In effect, root 
damage is lower when shearing process takes place. 

3. The engine speed has a stronger effect on the blade root damage especially at higher RPM (> 
5000). This is due to the fact centrifugal loading is ~ i2 2 and hence contributes more towards 
the root damage. 

4. The root damage increases linearly with the increase in impact radius along the blade and the 
slope of the line is very high at larger engine speeds. 
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Figure 1 . — A schematic of ice impact on an engine blade. 
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Figure 2.— Geometry of ice impact on an engine blade. 
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(a) Engine Speed * 3000 RPM 








Figure 6.— Effect of ice speed on blade root damage and Figure 7. — Effect of engine speed on blade root damage, 

impact angle. 



Engine Speed (RPM) ► Impact Location (% span) — ► 



Impact Angle, 0 O — ► 














REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response. Including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and completing and reviewing the collection of information Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services. Directorate tor information Operations and Reports, 1215 Jefferson 
Davis Highway. Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget. Paperwork Reduction Project (0704-01 B8), Washington, DC 20503. 


1. AGENCY USE ONLY ( Leave blank) 


4. TTTLE AND SUBTITLE 


2. REPORT DATE 


August 1992 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


Root Damage Analysis of Aircraft Engine Blade Subject to Ice Impact 


5. FUNDING NUMBERS 


6. AUTHOR (S) 

E.S. Reddy, G.H. Abumeri, C.C. Chamis, and P.L.N. Murthy 


WU -505-62-0K 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135-3191 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-7206 


8. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC. 20546-0001 


10. SPONSORING/MONrrORING 
AGENCY REPORT NUMBER 


NASA TM- 105779 


11. SUPPLEMENTARY NOTES 

E.S. Reddy and G.H. Abumeri, Sverdrup Technology, Inc., Lewis Research Center Group, 2001 Aerospace Parkway, 
Brook Park, Ohio 44142. C.C. Chamis and P.L.N. Murthy, NASA Lewis Research Center, Cleveland, Ohio. 
Responsible person, P.L.N. Murthy, (216) 433-3332. 


12s. DISTRIBUTION/ AVAILABILITY STATEMENT 


12b. DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 39 


13. ABSTRACT (Maximum 200 words) 

The blade root response due to ice impact on an engine blade is simulated using the NASA in-house code BLASIM. 
The ice piece is modeled as an equivalent spherical object impacting on the leading edge of the blade and has the 
velocity opposite to that of the aircraft with direction parallel to the engine axis. The effect of ice impact is consid- 
ered to be an impulse load on the blade with its amplitude computed based on the momentum transfer principle. The 
blade response due to the impact is carried out by modal superposition using the First three modes. The maximum 
dynamic stresses at the blade root are computed at the quarter cycle of the first natural frequency. A combined stress 
failure function based on modified distortion energy is used to study the spanwise bending damage response at the 
blade root. That damage function reaches maximum value for very low ice speeds and increases steeply with 
increases in engine speed. 


14. SUBJECT TERMS 


Engine blades; Ice impact modeling; Root damage analysis; Transient response 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 10. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


15. NUMBER OF PAGES 

12 


16. PRICE CODE 

A03 


20. LIMITATION OF ABSTRACT 


NSN 7540-01-280*5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std Z30-18 
208-102 


























